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Abstract 

We study the origin of the large thermopower in quasi-two-dimensional a r-type organic conductor, r- 
(EDO-S ) S-DMEDT-TTF) 2 (AuBr 2 ) 1+ , / (y < 0.875), from the view point of a "pudding mold" -type band 
structure. We calculate the electronic band structure using an ab initio band calculation package, and 
obtain a tight binding model fit to the ab initio band structure. Using the model and the Boltzmann's 
equation approach, we calculate the temperature dependence of the Seebeck coefficient. We conclude that 
the peculiar band structure is the origin of the large Seebeck coefficient and the appearance of the maximum 
value at a certain temperature. 

Keywords: r-type organic conductor, quasi-two-dimensional system, thermopower, band structure 



1. Introduction 

Recently, large thermopower has been exper- 
imentally observed in quasi-two-dimensional or- 
ganic conductors having r-type alignment of the 
donor molecules, e. g. r-(EDO-5",5"-DMEDT- 
TTF) 2 (AuBr 2 )i +2/ {y < 0.875)[l|. The absolute 
value of the Seebeck coefficient takes its maximum 
of about | S | ~ 150 /J.V/K around the temperature of 
T ~ 100 K to T ~ 150 K. 

A possible link between the large thermopower 
and the electronic structure is of special interest. 
In fact, a calculated band structure for r-(EDO- 
5 , ,5 , -DMEDT-TTF) 2 (AuBr 2 )i + y using the extended 
Hiickel method [3] has shown presence of a star- 
shaped Fermi surface and a peculiar band structure 
with flat portions near the T point. Another feature 
of this material is the variation of the band-filling 
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corresponding to the amount of anions. The anion 
content of y corresponds to the electron band filling 
of n = 1.5 - y/2. 

In the present study, we perform an ab initio band 
calculation for r-(EDO-5*,S'-DMEDT-TTF) 2 AuBr 2 
using the WIEN2K package. 5 and make a tight 
binding model fit to the obtained ab initio band 
structure. We calculate the Seebeck coefficient using 
the Boltzmann's equation approach and show that 
the peculiar band structure is the origin of the large 
Seebeck coefficient. 



2. Band structure 

We have performed calculations using all-electron 
full potential linearized augmented plane-wave 
(LAPW) + local orbitals (lo) method to solve the 
Kohn-Sham equations using density functional the- 
ory (DFT) within the framework of WIEN2K Q. 
This implements the DFT with different possible ap- 
proximation for the exchange correlation potentials. 
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The exchange correlation potential is calculated us- 
ing the generalized gradient approximation (GGA). 
The single-particle wave functions in the interstitial 
region are expanded by plane waves with a cutoff 
of i?MT-Kmax = 3.0, where Rmt denotes the small- 
est muffin tin radius and if max is the maximum 
value of the K vector in the plane wave expansion. 
In r-(EDO-£,£-DMEDT-TTF) 2 AuBr2, the muffin- 
tin radii are assumed to be 2.38, 2. If, 1.61, 1.27, 
1.18 and 0.64 atomic units (au) for Au, Br, S, O, 
C and H, respectively. For the value of K max and 
the plane wave cutoff energy, A max — 4.7 and the 
plane wave cutoff energy is 298.8 eV for r-(EDO- 
S,S-DMEDT-TTF) 2 AuBr 2 . The wave functions in 
the muffin tin spheres are expanded up to l m ax = 10, 
while the charge density are Fourier expanded up to 
G max = 20. Calculations are performed by using 512 
fc-points in the irreducible Brillouin zone. 

The calculated band structure is shown in Figure 
2] (a). We have made a tight binding model fit to the 
ab initio band as shown in Figure 0] (b) , where the 
transfer energies are determined as t\ = 161.61, t 2 = 
-157.29, t 3 = -11.44, £ 4 = 21.26, and t 5 = -2.81 
[meV]. As seen in Figure [4] (a), there are two band 
dispersions near the Fermi level, separated by a band 
gap. The band width is estimated to be 1.28 eV, and 
the maximum value of the band gap is around 0.f3 
eV. 

3. Boltzmann's equation approach 

Using the Boltzmann's equation approach, the 
thermopower is given by 

S = 1 K 1 , (f) 

where e(< 0) is the electron charge, T is the temper- 
ature, tensors Ko and Ki are given by 

K„ = £ r (fc) „ (k) v (k) ( \ {e k - M } w , 

k { k ) 

Here, e (k) is the band dispersion, v (k) = V^,£ (k) 
is the group velocity, r (fe) is the quasiparticle life- 
time, / (e) is the Fermi distribution function, and /i 



is the chemical potential. Hereafter, we simply re- 
fer to (K n ) xx as K n , and S xx = (l/eT)(Ki/K ) (for 
diagonal K ) as S. Using K , conductivity can be 
given as a xx = e 2 A'o = o. Roughly speaking for a 
constant r, 

K ~Y,\v 2 A + v%), K\ — (k B T)H'(v B - v\) (3) 

(apart from a constant factor) stand, where £' is a 
summation over the states in the range of \e(k) — 
fi\ <^ ksT, and va and vb are typical velocities for 
the states above and below li, respectively. In usual 
metals, where va ~ vb, the positive and negative 
contributions in K\ nearly cancel out to result in a 
small S. Now, let us consider a band that has a 
somewhat flat portion at the top (or the bottom), 
which sharply bends into a highly dispersive portion 
below (above). We will refer to this band structure 
as the "pudding mold" type [3| . For this type of band 
with \x sitting near the bending point, v\ 3> v B holds 
for high enough temperature, so that the cancellation 
in K\ is less effective, resulting in \K\\ ~ (fc^TjE'i)^ 
and Kq ~ ^'v 2 A , and thus large \S\ ~ 0(k B /\e\) ~ 
O(100)/iV/K. An important point for this type of 
band is that the large v A and the large Fermi surface 
result not only in large |^-| oc S but also in large 
K a oc a as well, being able to give a large power 
factor S 2 a. The pudding mold type band scenario 
well explains the coexistence of large thermopower 
and metallic conductivity observed in Na x Co02 

Now, if we look at the band structure of the r-type 
conductor from this viewpoint, the lower band has a 
nearly flat portion at the top, while the upper band 
has a nearly flat portion at the bottom. Therefore, 
the band structure can be considered as "inverted 
pudding mold type band" on top of a "pudding mold 
type band" . 

4. Thermopower 

(2) Using the Boltzmann's equation approach and as- 
suming a rigid band, we have calculated the tem- 
perature dependence of the Seebeck coefficient for 
the band filling n = 1.0625 and n = 1.125, which 
corresponds to the anion content of y — 0.875 and 
y = 0.75, respectively. Although the maximum ab- 
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Fig.l (a) The band structure for the ab initio result 
(red solid lines) and the tight binding fit (blue dashed 
curves), (b) The tight binding model in this study. 

solute value of the calculated S is somewhat smaller 
than the experimental one [l|, the overall tendency 
of the temperature dependence is well reproduced by 
the calculation when we adopt y — 0.875. 

The reason for the appearance of the maximum 
value of S can be understood as follows. In the low 
temperature regime, due to the presence of the band 
gap, only the upper band is effective, which gives a 
negative Seebeck coefficient. The Seebeck coefficient 
in this temperature regime can be large (despite the 
metallic conductivitiy) because of the inverted pud- 
ding mold shape of the upper band. At higher tem- 
peratures, the lower band, which gives a positive con- 
tribution to S, becomes effective, and this suppresses 



Fig. 2 The calculated Seebeck coefficient as functions 
of the temperature for the band filling of n — 1.0625 
(solid red) and 1.125 (dashed green), corresponding 
to y = 0.875 and 0.75, respectively. 

the absolute value, resulting in a maximum value at 
a certain temperature. 

5. Conclusion 

We have studied the origin of the large ther- 
mopower in the r-type organic conductor. We have 
performed an ab initio band calculation by using 
WIEN2K, and made a tight binding model fit to 
the ab initio band. Using this tight binding model, 
we have calculated the Seebeck coefficient which 
nicely reproduces the experimental results. We 
conclude that the peculiar band structure consisting 
of pudding + inverted pudding pudding mold type 
bands is the origin of the large thermopower as well 
as the appearance of the maximum value at a certain 
temperature. A more detailed analysis (such as on 
the molecule dependence) is now underway Q. 
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